Classical phase-space descriptions of continuous- variable teleportation 
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The nonnegative Wigner function of all quantum states involved in teleportation of Gaussian 
states using the standard continuous-variable teleportation protocol means that there is a local 
realistic phase-space description of the process. This includes the coherent states teleported up 
to now in experiments. We extend the phase-space description to teleportation of non-Gaussian 
states using the standard protocol and conclude that teleportation of non-Gaussian pure states with 
fidelity of 2/3 is a "gold standard" for this kind of teleportation. 



Quantum teleportation is a process that can trans- 
fer an arbitrary quantum state from a system held by 
one party, usually called Alice, to a system held by a 
second party, usually called Bob. The process requires 
a pair of systems, shared by Alice and Bob, in an en- 
tangled state — the entangled resource — and a "small" 
amount of classical information transmitted from Alice 
to Bob. Originally proposed for qubit states [1], telepor- 
tation protocols were later extended to states of a sys- 
tem described by continuous phase-space variables, such 
as a massive particle or a mode of an optical field [2, 3]. 
This continuous-variable teleportation protocol was im- 
plemented in an experiment that teleported a coherent 
state of an optical-frequency electromagnetic mode with 
fidelity 0.58 ± 0.02 [4]. Two recent experiments have im- 
proved the experimental fidelity of the teleported coher- 
ent state to values of 0.64 ± 0.02 [5] and 0.61 ± 0.02 [6]. 

In the standard continuous- variable teleportation pro- 
tocol [3], Alice and Bob share an entangled Gaussian 
state of two modes, A and B, which have annihilation 
operators a and b; this entangled resource is ideally a 
two- mode squeezed state [7]. The state to be teleported 
is the pure [8] state p — IV'XV'I of a mode V in Alice's pos- 
session, which has annihilation operator v. The protocol 
consists of (i) Alice's measuring two (commuting) joint 
quadrature components of modes V and A, specifically 
the Hermitian real and imaginary parts of the operator 
V + a', (ii) Alice's communicating the (complex) result 
£ to Bob, and (iii) Bob's displacing mode B by £. The 
efficacy of the protocol is quantified by the fidelity be- 
tween the output state of mode B and the input state 
\ip), averaged over the possible measurement results. 

Experiments to date have teleported only coherent 
states. It is generally believed, though not proved, that 
teleporting coherent states with average fidelity T > 1/2 
requires an entangled resource [9]. It has thus been ar- 
gued that teleportation of coherent states with fideli- 
ties above 1/2 constitutes truly quantum teleportation 
[9, 10]. Using a variety of arguments, other workers 
have contended that T = 2/3 is the appropriate bound- 
ary between classical and quantum teleportation [11, 12]. 
While acknowledging the need for an entangled resource 



for teleporting coherent states with fidelity T > 1/2 (in- 
deed, we provide additional evidence), we add a fresh 
perspective by investigating whether the entangled re- 
source is used in a way that can be accounted for by a 
classical phase-space description. When such a descrip- 
tion exists, it provides a local realistic hidden-variable 
model for the teleportation protocol. 

Our investigation is motivated by the fact that the co- 
herent states and the Gaussian entangled resource used 
in the experiments have nonnegative Wigner functions 
[13], which are phase-space probability distributions that 
provide a classical description of measurements of the 
quadrature components. A nonnegative Wigner func- 
tion does not give a classical description of measurements 
other than those of quadrature components; specific such 
measurements on an entangled Gaussian state cannot be 
given a local realistic description and thus violate Bell 
inequalities [14]. Since the standard protocol uses only 
quadrature measurements, however, we conclude that 
for teleporting coherent states — or any Gaussian state — 
using the standard teleportation protocol, the nonnegative 
Wigner function of the three modes gives a classical, lo- 
cal realistic description for all fidelities. This means that 
all the experiments to date — and any such experiment 
that teleports coherent states, no matter what fidelity is 
achieved — can be accounted for in terms of purely classi- 
cal correlations, with no need for a quantum-mechanical 
explanation. 

To find situations where the Wigner function does not 
provide a classical phase-space description of the stan- 
dard protocol, we must look to teleportation of non- 
Gaussian states, which (for pure states) have Wigner 
functions that take on negative values [15]. To accommo- 
date non-Gaussian states, we extend our hidden-variable 
model by allowing (i) Alice to substitute a randomly dis- 
placed state with nonnegative Wigner function in place 
of the non-Gaussian state and (ii) Alice and Bob to cheat 
by teleporting this new "smeared-out" state with perfect 
fidelity. Teleportation of non- Gaussian pure states with 
fidelity T > 2/3 cannot be accommodated within this ex- 
tended hidden-variable model, thus making a fidelity of 
2/3 a "gold standard" for teleportation of non-Gaussian 
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pure states. 

We begin with a brief Wigner-function-based review 
of the teleportation protocol. The state pab of modes 
A and B has Wigner function Wab{&, P), which is a 
quasidistribution for the c-number complex amplitudes 
a and (3 corresponding to the annihilation operators a 
and b [16]. In the standard protocol, Wab{&, (3) is a 
Gaussian, but for the present, we allow it to be a gen- 
eral Wigner function. The (pure [8]) state p = IV'XV'I 
of mode V has Wigner function W p {v), where v is the 
c-numbcr complex amplitude corresponding to annihila- 
tion operator v. The overall Wigner function of the three 
modes is W p {v)W AB (a, /?)■ The state of mode B after 
a measurement of v + a' that yields result £ has Wigner 
function 

W'(P\0 = ^J d 2 vd 2 a 5{u + a*- £) W p [v)W AB (a, 0) , 

(1) 

where 

p(0 = j d 2 V d 2 ad 2 p5(v + a*-£,)W p {v)W AB (a,l3) (2) 

is the probability to obtain result £. 

Having received result £ from Alice, Bob displaces 
the complex amplitude of mode B by yielding a 
state p O ut(0 with Wigner function W ou t(/?|£) = W'{j3 — 
The fidelity of this output state and the input 
state is F{£) = (V'|Pout(C)IV ; )- We are interested in 
the average of this fidelity over all measurement rc- 
suits, T = \d\y{i)F(t) = MPoutIV), where p out = 
J d 2 £p(£)p out ((;) is the average output state, having 
Wigner function 



W T>oM= d\p{0W o M = / d 2 vG{v)W p {(i - v) 



(3) 



Here 



G{v) = J d 2 ad 2 /38(l3 + a* -v)W AB (a,[3) (4) 

is the (nonnegative) probability to obtain result v in a 
measurement of b + a 1 * on modes A and B. Equation (3) 
shows that the average output state is a mixture of dis- 
placed input states, 



d 2 vG{v)D{v)pD\v) , 



(5) 



where D(v) is the displacement operator. 

We can now write the average output fidelity in two 
complementary forms, 

T = J d 2 vG{v)\C p {v)\ 2 

= it j ' d 2 pd 2 v G{l3~v)W p {!3)W p {v) , (6) 



where 

C p {v) = - J d 2 pW p (p)e^'-»'» , (7) 

the symmetrically ordered (Wigner- Weyl) characteristic 
function of the input state [17], is the Fourier transform 
of the Wigner function. The first form in Eq. (6) comes 
directly from Eq. (5), and the second from from writing 
the fidelity as an overlap of the Wigner functions for the 
input and average output states. The effect of the initial 
state of modes A and B on the average fidelity is con- 
tained wholly in the marginal distribution G(v). High- 
fidelity teleportation occurs when G{v) is very narrow, 
i.e., when the quadrature components contained in 6 + at 
are sharp, expressing a particular kind of correlation be- 
tween modes A and B. Using the Fourier transform (7), 
we can derive two other, equivalent forms for the average 
fidelity, 

T = tt J d 2 (3d 2 vG{p-v)W p {(3)W p {v) 
= J d 2 vG{ V )\C p { V )\ 2 , 

G{y) = J d 2 pG(p)e^'- v *» 



where 



(8) 



(9) 



is the Fourier transform of G(p). 

Before proceeding to the standard protocol and our 
hidden-variable models, we pause here to demonstrate 
the one technical result we need. We wish to find the 
maximum value of the integral 



I = J d 2 ad 2 pe- t ^-^ 2 / 2 W AB (a 1 P) , t > 0, 



(10) 



over the Wigner function Wab {pi, [3) of a joint state p AB 
of modes A and B. Introducing annihilation operators 
c = (a + b)/\/2 and d = {a — b)/V2, with corresponding 
c-number variables 7 and S, we can rewrite I as 

1=1 d 2 1 d 2 5e- t W 2 W CD {l,6) = J d^e-^WoiS) , 

(11) 

where Wcd{"I, 5) = Wab{oi,(3) is the Wigner function 
written in terms of modes C and D, and Wd{8) is the 
Wigner function for mode D alone. The integral can 
now be thought of as the expectation value, tv{A t po) = 
tv{A t pAB), of the D-mode operator A t whose symmetri- 
cally ordered associated function [17] is e - *' 5 ' , this op- 
erator being A t = {l + t/2)- 1 [{l - t/2)/(l + t/2)] dU . 

The integral / being the expectation value of A t , I is 
bounded above by the largest eigenvalue of A t . Since A t 
is diagonal in the number-state basis, with eigenvalues 
that decrease in magnitude with the number of quanta, 
we have 

1 



= (largest eigenvalue of A t ) = 



l + t/2 ' 



(12) 
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with the maximum achieved if and only if pab is the 
vacuum state for mode D. 

For the case that pab is a pure product state, |\&) = 
\"4>a) ® IV'b), which turns out to be the case of interest 
here, the condition for achieving the maximum becomes 
rfl*) = or, equivalently, o,\4>a) ® \iPb) = \i>A) ® b\ipB), 
from which it follows that a\ipA) — \ipA)(ipB\b\ipB) and 
b\ip B ) = \i)B){^A\a\i>A), implying that \tp A ) and \tp B ) 
are identical coherent states. Thus the only pure product 
states that achieve the maximum in Eq. (12) are products 
of identical coherent states. 

We can use Eq. (12) to get one interesting result im- 
mediately: The maximum average fidelity for tele-porting 
a coherent state using the standard protocol, but with a 
separable state for modes A and B, is 1/2. To show this, 
suppose first that modes A and B are in a pure product 
state, with factorizable Wigner function Wa(oi)Wb(/3). 
The characteristic function for any coherent state satis- 
fies |C co h(^)| 2 = e - !"! , so we can use Eq. (4) and the 
first form in Eq. (6) to write the average fidelity as 

T = J d 2 ad 2 /3 e-^-^ 2 W A (-a*)W B (l3) . (13) 

Here Wa{—(x*) is the Wigner function for the time- 
reversed, parity-inverted state of mode A. The t = 2 
general bound (12) implies T < 1/2, with equality if and 
only if mode A is in a coherent state |a) and mode B is in 
the time-reversed, parity-inverted coherent state |— a*). 

Now suppose modes A and B are initially in a sepa- 
rable state, thus having a pure product-state ensemble 
decomposition. The fidelity is the average over the pure 
product-state ensemble, which shows that the fidelity is 
still bounded above by 1/2, with equality if and only if 
the separable state is a mixture of product states of the 
form | a) ® |— a*). This does not show that 1/2 is the 
maximum fidelity for coherent-state teleportation in the 
absence of entanglement, since the result applies only to 
the standard protocol, but it is an additional piece of 
evidence, distinct from the results reported in Ref. [10]. 

We now take up again our analysis of the standard 
teleportation protocol, assuming that modes A and B 
are in a Gaussian state with Wigner function 



W AB (a,/3) 



4(c 2 - s 2 ) 



-2c(H 2 + |/?| 2 )-2s(a/3+ Q */3*) 



(14) 

where c and s satisfy |s| < c < Vl + s 2 . This state is 
pure if and only if c = \/l + s 2 , in which case the state 
becomes a two-mode squeezed state with c = cosh 2r and 
s = sinh2r, where r is the squeeze parameter [7]. The 
state (14) is separable if and only if c < 1 — |s| [18]. 

For the Wigner function (14), the distribution (4) is a 
Gaussian, 



W 7T 7rf 



where t = 2/(c + s) is the single parameter needed to 
characterize the fidelity that can be achieved with this 
entanglement resource. The Wigner function (3) of the 
average output state is the (s = — i)-ordered quasidistri- 
bution, Wp \v) [17], of the input state: 

W~ Po JP) = ^J d?ve-^-^l*W p {y) = W^(u) . 

(16) 

For t = 0, G{v) is a (5-function, and the output state 
is identical to the input state (perfect teleportation). 
For t = 1, the Wigner function of the average output 
state is the Husimi Q distribution of the input state, i.e., 

Wp mt (J3) = W^ifi) = Q„(J3) = (I3\p\P)/k- 

For < t < 2 (c > 1 - s), the state (14) is entangled, 
with the right sort of correlations for this protocol, these 
correlations decreasing as t increases. At t = 2, the state 
passes through the separability boundary c + s = 1 , and 
for t > 2, the state either is separable (c < 1 — |s|) or, 
though entangled (c > 1 + s), has the wrong sort of 
correlations for this protocol. 

The average fidelity of Eqs. (6) and (8) now becomes 



2 

Vt 



J dVe- 2 !"! 2 /*!^!/)! 2 
= - t J d 2 [3d 2 ue- 2 \l ) - v \ 2 ' t W p {r3)W p {u) 

= J d 2 (3d 2 V e-^-»\ 2 l 2 W p {l3)W P {v) 

= IJ d 2 ve-W 2 / 2 \C p (v)\ 2 . (17) 

The second form is the overlap of the Wigner function 
and the s-ordered quasidistribution for the input state. 
The first two derivatives of the last form show that T p (t) 
is a strictly decreasing, strictly concave function of t. 
These forms also show that the average fidelity obeys the 
scaling relation T p (t) — 2J r p (4 : /t)/t, which again draws 
attention to the separability boundary at t = 2. 

Given that |C co h(^)| 2 = e - '"' , the average fidelity for 
teleporting a coherent state is J r C oh(^) = (1 + i/2) -1 - 
For number states \n), whose Wigner functions take on 
negative values (except for n — 0), we have obtained an 
analytic formula for a generating function 



(18) 



n=0 



~ y/(l + t/2) 2 - 2A(1 + i 2 /4) + A 2 (l - t/2) 2 ' 
The resulting fidelity for teleporting a number state is 



(1 - t/2)" D /l+t 2 /4 



P, 



+ t/2) n + 1 n \l-t 2 /4 



(19) 



2M7* 



(15) 



where P n {x) is a Legendre polynomial. This 
gives F\ n)(n \{l) - 2P„(5/3)/3" +1 and F\ n){n \{2) = 
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(2n)!/2 2n+1 (n!) 2 . Another example of a state with a 
Wigner function that takes on negative values is the su- 
perposition \tp) = (|0) + |l))/-v/2, which has teleportation 
fidelity F ww (t) = (l + 3t/4 + t 2 /4)/(l + t/2) 3 . 

We now return to our Wigncr-function-based discus- 
sion of hidden-variable models for teleportation. For 
Gaussian input states and for the two-mode entangled re- 
source (14), all the Wigner functions are nonnegative [13], 
so they provide a classical phase-space description — and 
hence a local hidden-variable description — for this kind 
of teleportation, no matter what fidelity is achieved. The 
hidden variables are the quadrature components of all the 
modes, and the overall Wigner function is a probability 
distribution for these hidden variables. 

All non-Gaussian input pure states have Wigner func- 
tions that take on negative values (Hudson-Piquet the- 
orem [15]) and thus cannot be incorporated in the 
simple hidden-variable model. To see what can be 
achieved within a classical phase-space description, sup- 
pose that before performing the teleportation protocol, 
Alice "kicks" the input state p randomly in phase space. 
The random kick is described by a Gaussian so that the 
average state after the kick is 

(J = A f d^e-^'^D^pD^v) . (20) 

7rf J 

We choose the kicking strength t to be the minimum 
value necessary to make p' have a nonnegative Wigner 
function, thus giving p' a classical phase-space descrip- 
tion and allowing it to be incorporated it within our 
hidden-variable model. For all non-Gaussian states, this 
minimum kicking strength is one vacuum unit, i.e., t = 1 
[19], implying that p 1 is the state whose Wigner func- 
tion, W p '{v) = Wp^^iv) = Q P {v), is the Q function 
of the original state p. Further suppose that Alice and 
Bob cheat by teleporting p' with perfect fidelity. Then 
the fidelity of the overall process is the overlap of the 
Wigner and Q functions of p, i.e., the t = 1 fidelity (17) 
of the standard protocol. Notice that bigger kicks (t > 1) 
would give smaller fidelity, making clear why we choose 
the smallest kicking strength consistent with giving p' a 
nonnegative Wigner function. 

These considerations, coupled with wanting to know 
the maximum teleportation fidelity for a given entangled 
resource t, motivate us to find the maximum value of 
the average fidelity T p (t) over all input pure states p, 
t = 1 being the value relevant for our hidden-variable 
model. The task can be restated as finding the pure 
state p that maximizes the overlap of the Wigner func- 
tion and the ,s-ordcrcd quasidistribution Wp S \v). We 
apply the bound (12) to the third form of the average 
fidelity (17), in this case maximizing over pure prod- 
uct states p ® p. The resulting maximum is T max {t) = 
(1 + t/2) _1 , achieved if and only if p is a coherent state. 

Each non-Gaussian input state p has its own thresh- 
hold fidelity, T p (l) < ^ r ma x(l) = 2/3, below which its 



teleportation can be accommodated within our extended 
phase-space hidden-variable model and above which it 
cannot. Thus teleportation of a non-Gaussian state with 
fidelity exceeding T P {1) is required to rule out an expla- 
nation in terms of classical phase-space correlations. A fi- 
delity of 2/3 emerges as a "gold standard" for continuous- 
variable teleportation in the sense that teleportation of 
any non-Gaussian pure state with T > 2/3 cannot be 
fitted within our extended hidden-variable model. This 
conclusion applies only to our phase-space-based hidden- 
variable model; we have not shown that there is no local 
hidden-variable model that can accommodate teleporta- 
tion fidelities of 2/3 or above. 
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